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Recently I proposed a simple dynamical network model for discrete space-time which self-organizes 
as a graph with Hausdorff dimension dn = 4. The model has a geometric quantum phase transition 
with disorder parameter {dn — dg) where dg is the spectral dimension of the dynamical graph. Self¬ 
organization in this network model is based on a competition between a ferromagnetic Ising model 
for vertices and an antiferromagnetic Ising model for edges. In this paper I solve a toy version of 
this model defined on a bipartite graph in the mean field approximation. I show that the geometric 
phase transition corresponds exactly to the antiferromagnetic transition for edges, the dimensional 
disorder parameter of the former being mapped to the staggered magnetization order parameter 
of the latter. The model has a critical point with long-range correlations between edges, where a 
continuum random geometry can be defined, exactly as in Kazakov’s famed 2D random lattice Ising 
model but now in any number of dimensions. 

PACS numbers: 04.60.-m; 89.75-Hc 


I. INTRODUCTION 

Networks play a role in several of the major contenders 
for a theory of quantum geometries. In the spin foam ap¬ 
proach [I| to loop quantum gravity Q spin networks rep¬ 
resent the quantum states of the gravity on a given man¬ 
ifold. In the causal dynamical triangulation approach 
(CDT) Q, triangulations of space-time play the role of a 
lattice regularization of the Einstein-Hilbert action used 
to search for a non-trivial ultraviolet (UV) renormal¬ 
ization group fixed point defining quantum gravity in a 
non-perturbative sense, an Mproach that goes under the 
name of asymptotic safety 0. Causal sets Q and ener¬ 
getic causal sets @ are attempt to build geometries from 
networks of causal relations. Quantum graphity mod¬ 
els Q attempt the same starting from graphs. Finally, a 
very recent approach Q is the idea of considering discrete 
space-time as a quantum network, obtained by growing 
simplicial complexes. 

In a recent publication Q I proposed a different ap¬ 
proach to quantum gravity in which networks play an 
even greater role. The idea is to turn things around 
and, instead of starting from the correct infrared (IR) 
geometric variables and use discrete space-time just as a 
regularization artefact, as in the CDT approach, to posit 
purely combinatoric information bits as the fundamental 
UV quantum variables of gravity and have these self- 
organize so that geometry and general relativity emerge. 

Of course this presupposes the existence of a model 
having an UV fixed point corresponding to disordered 
bits and an IR limit in which vertices and edges self- 
organize to form a network with the topology of dis¬ 
crete space-time. In @ I have proposed exactly such 
a model, suitably generalizing Kazakov’s famed random 
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lattice Ising model 0. As I have shown in Q there is 
ample numerical and analytical evidence supporting the 
fact that the Hausdorff dimension of the emergent space- 
time in this model is predicted to be 4, the upper critical 
dimension of the Ising model [HI- 

This model distinguishes itself from previous ap¬ 
proaches in two major respects. First of all, it is much 
more radical, in the sense that no constraints at all, like 
the requirements of triangulations, causal sprinklings of 
Lorentz manifolds or simplicial complexes, are imposed 
on the discrete basic structure: this is completely deter¬ 
mined by self-organization driven by a minimum energy 
principle alone. Secondly, nothing is assumed a priori 
about time: contrary to most previous authors I believe 
time should not be assumed by imposing causality from 
the very beginning, but should be, rather, ’’explained” 
by the model itself. Admittedly, though, I am still quite 
short of this goal. The idea, however, is to explore if 
space-time can emerge from the simplest combinatoric 
variables with no other driving force than the competi¬ 
tion arising from frustrated Ising Hamiltonians. 

II. THE MODEL 

The model is formulated in terms of N spin 1/2 bits 
Si = ± 1 , i = whose values indicate the pres¬ 

ence/absence of (Euclidean) space-time, respectively and 
N{N — l)/2 symmetric spin 1/2 bits Wij = wji = 0,1, 
i,j = 1... N, whose values denote a connection (1) be¬ 
tween space-time vertices Si and Sj or the absence thereof 

( 0 ) {wii = 0 ) 

'^ik'^kj 2 ( 1 ) 

S’ 

where J is the unique dimensionless coupling. I use units 
in which h = 1, c = 1 and all energies are measured in 
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units of the standard Ising coupling, second term in ([T]), 
which is set to one for simplicity of presentation. 

The second term in this energy function is the standard 
ferromagnetic Ising model. If the links Wij would be uni¬ 
formly drawn from random adjacency matrices of degree 
4, the model would be exactly Kazakov’s random lattice 
Ising model in two dimensions M- The first term in the 
energy function, on the other side, is simply a nearest- 
neighbours (sharing a common vertex) antiferromagnetic 
Ising model for the link spins. The generalizations with 
respect to Kazakov’s model, thus consist in dropping the 
restriction to degree 4 and drawing the random adjacency 
matrices from a Gaussian distribution. 

The driver of self-organization in this model is the 
competition between the vertex ferromagnetic coupling 
and the link antiferromagnetic one, creating ’’link frustra¬ 
tion” . Indeed the vertex ferromagnetic coupling favours 
the creation of many links (in a vertex aligned configu¬ 
ration) while the antiferromagnetic link coupling tends 
to suppress such links. As I have shown in [9|, for 
J = I/(4d— I), with integer d, the compromise is a 
2fi-regular ground state graph with power-law extension, 
exactly what one would expect for a discretized space- 
time. The spectral and Hausdorff dimensions [l^ of this 
“space-time graph” are determined by the unique dimen¬ 
sionless coupling J of the model. The spectral dimension 
is simply d while all evidence supports the fixed value 
4 for the Hausdorff dimension dn in the whole range of 
couplings between the two quantum phase transitions at 
d = 1, corresponding to the lower critical dimension of 
the Ising model, where space-time vertices themselves be¬ 
come disordered and d = 4, the upper critical dimension 
of the Ising model where spectral and Hausdorff dimen¬ 
sion start to coincide. 


III. A TOY VERSION OF THE MODEL ON 
BIPARTITE GRAPHS 

The purpose of the present paper is to solve a sim¬ 
plified version of © in the mean field approximation to 
elucidate the nature of the upper phase transition with 
disorder parameter (dn — d) and to show how critical 
graphs with long-range correlations and power-law ex¬ 
tension emerge as a result of frustrated edge antiferro¬ 
magnetism. 

Near the upper transition all space-time vertices are 
aligned. I will thus assume a configuration with Si = -1-1, 
VL In this case the model o reduces to a dynamical 
graph problem with an edge Hamiltonian 

z !:<'.(.). (2) 

» e(z)//(i) i e(i) 

where e{i) = 0,1 and f{i) = 0,1 denote edges emanat¬ 
ing from vertex i. In this case the frustration reduces 
to the competition between the edge antiferromagnetic 
interaction (first term) and an external edge “magnetic 


field” generated by the aligned space-time vertices (sec¬ 
ond term). In this space-time-aligned phase, the model 
can be thought of as being defined on the complete graph 
on N vertices, the (-1-1) edges defining a dynamically gen¬ 
erated subgraph representing the emergent space-time. 
The connectivity of this subgraph is dynamically deter¬ 
mined by the coupling constant J: for J = l/(4d— 1) 
with integer d the subgraph is 2d-regular. The toy model 
I will consider in this paper is defined by the same Hamil¬ 
tonian ([2]) but restricted on a connected, bipartite 2d- 
regular graph, while J is left free. As I will now show, 
this model is solvable in the mean field approximation. 

As a first step I will introduce a more familiar notation 
by defining standard spins ^e(i) = ±1 as ^e(z) = ~ 1- 

The Hamiltonian ([2]) reduces then to 

^ “ g"E E ~ ^EE^®(d ’ 

k (3) 

4 

The second step is to invoke Kdnig’s edge colouring 
theorem [l^ , which states that the edge chromatic num¬ 
ber of any bipartite graph equals its maximum vertex 
degree, in this case 2d since the graph is assumed 2d- 
regular. This means that at every vertex one can colour 
the incident edges with exactly 2d different colours, with¬ 
out two edges of the same colour ever touching. Let me 
now further subdivide the 2d colours into two types of 
light and dark colours: at every vertex there will be ex¬ 
actly d light colours and d dark colours. Exactly like one 
standardly treats antiferromagnets on bipartite lattices 
by introducing two different magnetizations for the two 
sub-lattices I will consider configurations with two differ¬ 
ent magnetizations mi and (defined as usual between 
-1 and 1) for the two types of edge colours and define 


m = 


md + mi 
2 


ms 



( 4 ) 


as the magnetization and the staggered magnetization, 
respectively. The original 2(i-regular graph corresponds 
clearly to a maximally “ferromagnetic configuration” 
TO = 1, TOs = 0. The maximally “antiferromagnetic con¬ 
figuration” TO = 0, TOs = 1, instead corresponds to a reg¬ 
ular connected sub-graph. In between these two extremes 
there are intermediate configurations in which each dark- 
coloured edge has a probability pd = {I + mi){2 of being 
present in the graph, whereas the probability for light- 
coloured edges is p; = (I -I- TOi)/2. 


IV. SOLUTION OF THE TOY MODEL IN THE 
MEAN FIELD APPROXIMATION 

In the mean field approximation I will decompose each 
edge spin into its mean value TOe(i) and a fluctuation 
around it: is{i) = We(i) -b(5^e(i), with = Cep) 
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and neglect terms quadratic in the fluctuations in the gives then two coupled equations for the two order pa- 
Hamiltonian, so that rameters m and rus in which one of them has no explicit 

dependency on the magnetic held, 


P'Hmf 


® e(i)^/(z) 

( 5 ) 

i e{i) 


2ms 

1 + TOg — 

2m 

1 + — TOg 


= th[4J'ms] , 

= th [2H - 4:J{2d - l)m] 


( 12 ) 


where /3 = 1/T is the inverse temperature and 

J = /3^ , H = ph. (6) 

At this point one can compute exactly the free energy 
per edge (a 2c?-regular graph on N vertices has dN edges 
by the degree sum formula). 

= -^ E • ( 7 ) 

{e} 

I will consider first the Ansatz of configurations with ex¬ 
actly d edges of mean and d edges of mean mi at 
every vertex. Apart from an irrelevant constant, in this 
case the free energy becomes 

fMF =-{d-l){ml + m'f)-2d mami 

— En cosh [2J {{d — l)md + dmi) — H] 

\J 

-En cosh [2J {{d - l)mi + dmd) — H] , (8) 

sj 

Note that, contrary to the usual situation in staggered 
antiferromagnets, where only vertex-spins of the two dif¬ 
ferent sub-lattices interact, here there are both interac¬ 
tions between the two types of light- and dark-coloured 
edges and amongst each type of colour also. Using O 
the free energy can also be easily rewritten in terms of the 
magnetization m and the staggered magnetization m^, 

Jmf =—{4:d — 2)m'^ + 2ml 

-In cosh [2((2d — l)m -|- mg) — H] 

— En cosh ((2d — l)m; — mg) — H] , (9) 

u 

The stationarity conditions dfMF/dm = 0 and 
dfMF/dmg = 0 lead to the equations 

TO-I-TOs = th [id — 2J7 ((2d — 1 )to — TOs)] , 
m — mg = th [id — 2jd ((2d — 1 )to-I-TO g)] . (10) 


Using the following formulas for the hyperbolic tangent 


th(x -d y) 
th(a; - y) 


th(a;) + th{y) 

1-d th(a:)th(2/) ’ 
th(a;) - th{y) 

1 — th(x)th(y) 


( 11 ) 


The first immediate observation about these equations is 
that, for id ^ 0 and finite Jd, the solution of the second 
equation always implies to ^ 0. Let me now consider 
the first equation for the staggered magnetization. The 
function on the left-hand side has a slope 2/(1 — mf ) at 
the origin, has two extrema at ±\/l — rrP and vanishes 
at ± 00 . The hyperbolic tangent on the right-hand side, 
instead has slope 4,J at the origin and approaches ±1 
at infinity. The important point is that the values of the 
left-hand side function at mg = ±1 are ±2/(2 —to^) which 
lie above (respectively below for negative mg) the hyper¬ 
bolic tangent asymptote. This means that, for small beta 
satisfying AJ < 2/(1 — m?) the only solution to the first 
equation is mg = 0. For J > 1/[2(1 — m?)], instead, two 
new solutions mg 0 appear. As always, one expects 
the solution mg = 0 to become unstable at this point. 

To check this point let me compute the second-order 
partial derivatives of the free energy. 


d^fMF 

dm^ 

d'^ fMF 

dmi 

d^fMF 

dmdmg 


= -4(2d - 1) - AJ{2d - 1)2 {1 - th2(-)} 

-4:d(2d-l)2{l-th2(+)} , 

= 4-4^{l-th2(-)}-4^{l-th2(+)} , 

= 4:d(2d-l){l-th2(-)} 
-AJ{2d-l){l-th^{+)] . (13) 


where I have introduced the short-hand notation 


th(±) = th [2jd ((2d — 1 )to ± TOg) — dd[ , 
th(—) = th [2jd ((2d — 1 )to — TOg) — dd[ . (14) 

Since d'^ fMF /dmf < 0, the sign of the determinant 
{d‘^fMF/dm‘^){d‘^fMF/dm1) - {d"^fMF/dmdmg)'^ is de¬ 
termined entirely by d^fMF / dmi.. Using m we can 
rewrite this as 


d"^ fMF 

dmi 


{mg 


0) = 4 - 8jd (1 - to2) 


(15) 


This shows that on the line ff = 1/[2(1 — to^)], 
where the two new solutions mg ^ 0 appear, 

the term d"^fMF/dmi. changes sign from positive to 
negative when ff is increased. Since, moreover 
{d^fMF/dmdmg) {mg = 0) = 0, it is the whole determi¬ 
nant {d'^fMF/dmf){d'^fMF/dm1) - {di^fMF/dmdmg)'^ 
that changes sign from negative to positive on the 
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line J = 1/[2(1 — m?)] when is increased. Fi¬ 
nally, since both {d‘^fMF/drn^){ms = 0) < 0 and 

{d'^fMF/drri^) {ms = 0) < 0 for J” > 1/[2(1—m^)] all this 
shows that, indeed, the unique high-temperature solution 
TOs = 0 becomes a local maximum when the two new so¬ 
lutions ms ^ 0 appear: these are thus the new stable so¬ 
lutions and at the critical temperature J7c = 1/[2(1 —m^)] 
there is a phase transition from a high-temperature phase 
characterised by to ^ 0 and TOs = 0 to a low-temperature 
phase with to ^ 0 and rUg ^ 0. 

We have already seen (eq. (fT31) ') that the second deriva¬ 
tive of the free energy vanishes at the transition. Let me 
further compute the third- and fourth-order derivatives: 


h 



d^fMF ^ _ th2( + )) - th(-) (1 - th^{- 

= 32J3 [(1 _ th2(+))^ + (1 - th2(-))' 


dml 

d^fmF 


dm^ 

+64J3 [th"(-) (1 - th"(-)) - th2(+) (1 - th2(+))] . 
Using again eq. m one can easily conclude that 
Jmf 


dml 

dmi 


{ms = 0 ) =0 


(TOg=0) =32J^ {1-m^f >0 , (17) 


which shows that, for finite J' (for which also to < 1) 
the phase transition is of second-order, i.e. continuous, 
although one cannot exclude that this is an artefact of 
the mean field approximation. 

The phase boundary can be derived easily by using 
eqs. m for TOg = 0 together with Jc = 1/[2(1 — m^)] 
and the inverse of the hyperbolic tangent: th“^(x) = 
(l/2)ln((l + l)/(l-x)), 


1 1 _|_ 'fll 

H = (4(i — 2)Jm H—In- 

2 1 — TO 


TO = a/1 - 


1 


or, reintroducing the coupling constant J of the original 
model. 


2d — \ T \ m 

h = - Jm H-In- , 

4 2 1-TO ’ 


4T 

w = \/l - — . (19) 


FIG. 1: The phase boundary of the antiferromagnetic transi- 
■)jjioin in the h-T plane, together with the fixed magnetic field 
of the model, eq. Q, for d = 4 and J = 1/10. The values 
hi {i = 1... 4) represent the critical magnetic fields where i 
edges decouple for each vertex at T = 0. 

(16) 

When the magnetic field is smaller than a second crit¬ 
ical value /i 4 = J/4, instead, the mean field equations 
OT imply a solution to = 0, TOg = 1. In the language of 
dynamical graphs this indicates a ground state consisting 
of a d-regular sub-graph that touches every vertex of the 
original 2d-regular graph on which the model is defined. 
This can be viewed either as a higher-dimensional gener¬ 


alization of a Hamiltonian cycle [1^ or as a self-avoiding 
hypersurface on the graph [l5l| . In the simple example of 
a hypercubic lattice, the spectral dimension is now clearly 
decreased to d/2, since a random walker has maximally d 
edges to choose from at every step, dg = d/2. The Haus- 
dorff dimension of this dynamical sub-graph, instead re¬ 
mains dH = d since the sub-graph visits every vertex of 
the original graph (exactly once), i.e. it is space-filling in 
the original graph. 

Between the two critical values = J/4 and he = 
{2d — l)T/4 the mean field equations (ITUl) have no solu¬ 
tion at r = 0. This indicates the failure, in this region, of 
the Ansatz of a configuration with half the edges of one 
type and half of the other type. It is easy, though, to re¬ 
peat the mean field computation with a generic Ansatz in 
which, at every vertex, there are na dark-coloured edges 
and ni light-coloured edges, with Ud + ni = 2d. The mean 
field equations become 


The phase boundary (fTOl) in the h-T plane, together 
with the magnetic field value implied by eq. m , are 
shown in Fig. 1 for d = 4 and J = 1/10. 

At zero temperature and above the critical magnetic 
field he = {2d — l)T/4 (= 0.175 for d = 4 and J = 1/10), 
the mean field solution implies to = 1 and Wg = 0. As 
already stressed above, this means that the ground-state 
configuration is the entire graph on which the model is 
defined. To make contact with geometry, let me consider 
the simple example of a d-dimensional hypercubic lattice. 
In this case, the spectral dimension dg coincides with the 
Hausdorff dimension dn, both being determined by the 
graph connectivity, dg = dn = d. 


TO + TOg = th [id — 2J {{2d — 1)to — (1 — {ud — n/))TOg)] , 
TO — ms = th [id — 2J7 ((2d — 1 )to -I- (1 -I- {ud — ni))ms(^0) 

While the upper critical magnetic field he = {2d — l)J/4 
remains unchanged, the lower value is now increased to 
hi = {1 {2d — 2i))J/4 for Ud = 2d — i and ni = i, i = 

I... d. This coincides with the upper value (2d — 1) J/4 
for ni = i = 1, rid = 2d— I : hi = he. This shows that, at 
T — 0, the quantum phase transition is of first-order, im¬ 
plying a jump from to = I, TOg = 0 to to = 0, TOg = I and 
the ’’decoupling of one edge”, since n; = I and to = 0, 

TOg = I together imply that at each vertex one edge of a 
given (light) colour is absent while all others are present. 
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In the hypercubic lattice example, this defines a regular 
sub-graph of spectral dimension dg = d — \/2 and Haus- 
dorff dimension dn = d. It is easy to show that, low¬ 
ering the magnetic field below the further critical value 
/i 2 = {2d — 3) J/4 at T = 0, the free energy is minimized 
by the Ansatz Ud = 2d — 2, m = 2. At this value an en¬ 
tire dimension (two edges) decouples. When h is further 
decreased, at T = 0, more dimensions decouple at hi, 
i = 3... d, until, at /i = J/4 the spectral dimension hits 
its lowest value d/2, as described above. It is to be ex¬ 
pected that these first-order quantum phase transitions 
extend as phase boundaries in the T > 0 region all the 
way to the mg —>■ 0 phase line. Unfortunately, the anal¬ 
ysis of these phase boundaries for T > 0 is much more 
complex and goes beyond the scope of the present paper. 
The important point here is, rather, that this toy model 
clearly shows how the statistical antiferromagnetic order 
parameter iris is related to the geometric disorder param¬ 
eter (dn — ds): when mg vanishes du — dg also vanishes, 
conversely, a non-zero value of rris indicates a decoupling 
between dg and dn- 

For the values of d and J shown in Fig.I, the transi¬ 
tion happens at finite temperature, since the magnetic 
field implied by the original model is a not a free 
variable but it is, rather, fixed itself in terms of the cou¬ 
pling constant J as in ([3]). This value is shown as an 
horizontal line in Fig. 1. The geometric transition de¬ 
scribed by the antiferromagnetic order parameter rUg is 
thus a second-order (in mean field theory) transition at 
a critical finite temperature Tg- At this temperature the 
dynamical graph becomes critical, with long-range corre¬ 
lations between edges on the distance measure defined by 
the embedding graph distance. The idea of the original 
model is that such a critical point defines a continuum 
space-time, exactly as in the original Kazakov model 0 
but not anymore restricted to 2D. 

Inserting the magnetic field ([3]) into eq. dD one can 
derive the phase boundary in terms of the original cou¬ 
pling constant J, 

J = , , I -r(l-2ri„i±I^) , 

(2c? — 1)(1 -I- m) \ 1 — mJ 

I 4T 

m =yi-— • (21) 

This shows that, for J = l/(4d — 2), the critical tem¬ 


perature becomes T = 0 and the geometric critical point 
corresponds to a purely quantum phase transition, as is 
shown in Fig. 2. This is exactly as in the full model Q 
with the only difference that there J = l/{Ad— 1) and 
the quantum phase transition corresponds to J = 4. In 
the present toy model this quantum phase transition is 
of first order, as explained above. Even a small tempera¬ 
ture, however, is sufficient to smoothen out the transition 
and make it of second-order (in the mean field approxi¬ 
mation) . 

Of course, the original model Q is defined on the 
complete graph K^, to allow every possible sub-graph 
h 



FIG. 2: The phase boundary of the antiferromagnetic transi¬ 
tion in the h-T plane, together with the fixed magnetic field 
of the model, eq. ((Sj, for d — 4 and J = l/(4d — 2) = 1/14. 
For this value of the coupling constant the geometric critical 
point corresponds to a purely quantum phase transition. 


as a dynamical ’’space-time graph” and, unfortunately, 
the complete graph is not bipartite. As in usual anti- 
ferromagnets, the frustration relations on non-bipartite 
graphs are typically much subtler than on bipartite 
graphs and lead to much more complex behaviour. It 
is possible that the recently developed network tensor 
methods [l3] may shed some light on the critical be¬ 
haviour of the fully frustrated model. For the moment, 
though, only numerical evidence is available in this case 
Q. This numerical evidence, however, supported by re¬ 
cent analytical results on finite-size scaling in the Ising 
model [ig points indeed to the existence of a geometric 
quantum phase transition at dn = Js = 4. 
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